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Abstract.  In  this  paper  a  new  transformation  referred  to  as  the  G^ ^-transform  is  introduced  for 

the  purpose  of  increasing  the  rate  of  convergence  of  a  sequence  to  its  limit.  In  particular  if  the 
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sequence  is  a  sequence  of  partial  integrals,  F(x),  then  G^ n '[F(x)]  is  shown  to  converge  super  fast  to 

(in} 

F(oo)  under  very  general  conditions.  The  G'  n  -transform  is  shown  to  be  closely  related  to  Levin  and 
Sidi’s  D-transform  and  the  Bn-transform  introduced  by  Gray,  et  al.  Several  examples  are  given. 
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1.  Introduction 


In  [2]  Gray,  Atchinson  and  McWilliams  introduced  the  Gm-transform.  They  showed  the 
transform  to  be  of  value  in  evaluating  improper  integrals,  especially  in  the  case  of  integrals  in  which 
the  integrand  decayed  like  an  exponential  function. 

The  results  of  Gray,  et  al  were  significantly  extended  by  Levin  and  Sidi  [6]  to  a  much  more 
general  case,  which  no  longer  required  the  exponential  behavior  of  the  integrand  in  the  limit  to  be 
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effective.  This  was  accomplished  by  introducing  what  they  referred  to  as  the  D-transformation.  Given 
a  sequence  of  values  of  a  function,  the  D-transformation  is  a  very  effective  tool  for  increasing  the  rate 
of  convergence  of  the  sequence  to  its  limit  in  most  cases.  In  fact,  to  date,  the  D-transformation  (or  its 
discrete  counterpart,  the  d-transformation)  is  probably  the  most  effective  general  method  for  increasing 
the  rate  of  convergence  of  a  sequence  available.  For  some  cases  its  effectiveness  has  been  further 
improved  recently  by  a  slight  extension  by  Gray  and  Wang  [4]. 

Since  the  original  paper  by  Levin  and  Sidi  a  number  of  papers  have  been  written  on  the  Levin 
and  Sidi  transformation  and  Ford  and  Sidi  [1]  have  shown  that  the  generalized  Richardson 
extrapolation  procedure  (GREP)  can  be  used  to  calculate  the  D-  and  d-transformations.  One  aspect  of 
the  D-transformation  that  has  not  been  investigated  very  closely  is  the  behavior  of  the  transformation 
as  the  step  size  between  an  equally  spaced  sequence  of  points,  where  the  functional  values  are  taken, 
approaches  zero.  In  the  particular  case  where  the  sequence  is  a  sequence  of  improper  integral 
approximations  and  the  integrand  is  given  at  equally  spaced  points,  i.e.,  at  Xq  4-  kh,  h  fixed  and  k  =  0, 
I,  2,  .  .  .  ,  N,  this  phenomonon  was  investigated,  [2],  for  the  Gm-transform.  In  particular  it  was 
demonstrated  that  in  the  case  of  the  Gm-transformation,  it  is  often  the  case  that  the  approximation 
improved  as  h— *0.  Investigating  this  property  then  led  to  a  new  transformation  referred  to  as  the  Bin- 
transformation  [3].  To  be  specific  let 


F(x) 


-  f  <w 


du 


(i) 


where  F(x)  — ►  S  <  oo  as  x  — *  oo.  Then  lim  ,  rF(x+kh)1  —  Bm[(F(x)].  The  Bm-transformation 

h— 0  J 

was  then  shown  to  be  of  value  in  its  own  right,  particularly  in  the  area  of  producing  functional 
approximations  to  tail  probabilities  (see  [3]).  The  fundamental  difference  between  Gm  and  Bin- 
transformations  is  that  the  former  requires  values  of  the  function  being  transformed  at  several  points 
while  the  Bm  requires  values  of  the  function  and  its  derivatives  at  a  single  point.  In  the  case  of  tail 


2 


probabilities  this  distinction  allows  one  to  replace  the  problem  of  integration  completely  by 
differentiation.  ,  , 

y  In  this  paper  we  introduce  a  new  transformation,  referred  to  as  the  transformation,  which  is 

in  fact  the  limit  of  the  Levin  and  Sidi  D- transformation  when  applied  to  a  function  evaluated  at 

/  \ 

equally  spaced  points.  Thus  the  Gn  -transformation  can  be  viewed  as  the  extension  of  the  Bm- 

transformation  analogous  to  the  Levin-Sidi  generalization  of  the  Gm-transformation.  In  fact,  as  we 

shall  see,  Bm  =  Gj  .  'The  Levin-Sidi  generalization  has  been  clearly  shown  to  be  a  powerful  tool  for 

evaluating  improper  integrals  when  the  integrand  is  given  at  a  sequence  of  pointsf  it  should  be 

(m) 

expected,  therefore  that  the  Gq  -transformation  will  be  useful  when  the  integrand  and  its  derivations 
are  given  at  a  single  point.  This  is  shown  to  be  the  cae  in  boththeory'  and  practice.  In  particular, 
under  very  general  conditions,  Gnm^[F(x)]  is  shown  to  converge  to  S  faster  than  F(x)  at  a  rate 

,  t 

proportional  to  x"n.  Several  examples  are  given.  - —  -  /  .  J,r  ,  - 


2.  Background  and  Definitions 
Let  F(x)  and  S  be  defined  as  in  equation  (1)  and  let 


X  •  *1 


fOO 

e(x)  =  S  —  F(x)  =  /  f(u)  du. 


Further  let 


UJx) 


l\,  oo 


=  x  £  Qk  i  x 

i=l  K,‘ 


— i+1 


(2) 


(3) 


where  o^  j  ^  0  and  is  an  integer  with  <  k.  Also  suppose  m  is  the  smallest  integer  such  that 
e(x)  satisfies  the  differential  equation 


Um(x)  y(m)  +  Um_1(x)  y(m  1}  +  .  .  .  +  Uj(x)  y'-y  =  0, 


(4) 


for  some  set  of  U^'s.  For  motivational  purposes  suppose 
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a,  •  =  Or0  ■  .  .  .  =  a  •  =  0 
1,1  2,1  ITU 


when  i  >  n.  Then 


i=l 


x*m-i+l  f(m) 


(x)  +  .  .  .  +  £  ai . 
i=l 


a.  €,  — i-t-1  , 

x  1  e^x)  — f(x)  =  0  . 


Thus 


F(x)  =  S  +  £  ,  x*1  i+1  f(x)  +  .  .  .  +  £  «ra .  x;”-1+1  tm~%) 

i=l  ’  i=l 


(5) 


and 


/k)(x)  =  ±  {x‘>-i+1  „P  +  . . .  +  £  Vi  { 


(k+1) 


(6) 


k  =  0,  1,  .  .  .  ,  N— 1,  N  =  mn  .  Now  let 


Hk+l(zj  5  aij)  = 


Z1  z2  *  • '  zk+l 


all  a12  '  •  •  al,k+l 


ak,l  ak,2  •  •  •  *k,k+l 


(7) 


and  define 


;[fw, 


G  F(x),  f(x),  .  .  .  ,  fv  '(x) 


I  aij(x)J  = 


HN+l(cj  ;  aij) 


(8) 


where  =  1,  C2  =  c-j  =  .  .  .  =  =  0  ,  z-  =  F^  *\x),  and  for  j  =  1,  2,  .  .  .  ,  mn+1, 
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£,-i+l 


iG-i) 


/  r  ^2~  i+°+i .  ~i 

ai  ;(x)  =  ( L*  2  f'wj 


G-i) 


£m-i+(m-l)n+l  itn- 


a-i) 


i  =  1 . n, 


i  =  n+1,  .  .  .  ,  2n, 


i  =  (m-l)n+l,  .  .  .  ,  mn 


From  equations  (5)  and  (6)  it  follows  from  Kramer’s  rule,  if  ;  a- j)  ^  0,  that 

g[f(x),  f(x),  .  .  .  ,  f(N_1)(x);  aij(x)]  =  S  =  J°°  f(u)du  . 


Now  suppose  f(x)  satisfies  (4),  i.e.. 


f(x)  =  £  ukw  f(k)w  > 

k=l 


where  Ujc(x)  is  defined  by  (3)  and  we  no  longer  assume  that  the  series  in  (3)  is  of  finite  length  but 
rather  that  it  converges  for  all  x  >  a.  Suppose  further  that 


xi^oo  uk^  ^  W  M  =  °> 


for  k=i,  t-4-1,  .  .  .  m,  i=l,  .  .  .  ,  m.  Then  Levin  and  Sidi  [6]  have  shown  that  there  exists  a  set  of  af  • 

K,l 


such  that 


umW  ^(x)  +  U^_j(x)  €^m  *\x)  +  .  .  .  +  U|(x)  e'(x)-e(x)  =  rn(x)  , 


where 


5 


(13) 


ukW  =  £  °k,i x 

i=l 


£k_i+1 


and 


r„(x)  =  £  f‘k'1)(x)  £  .J4  x'k  i+1  =  o(x-"  t(x))  , 


(14) 


k=l 


i=n+l 


m  e.  /k  — 1) 

since  t(x)  -  £  a?  ,  x  F  (x).  The  notation  implies  that  if  p(x)  -  q(x)  then  p(x)/q(x)  — ►  1, 
k=l  ’ 


as  x  — *  oo  .  From  (12)  we  have 


F(x)  =  S  +  £  Uj(x)  f(k  1}(x)  +  rn(x), 
k=l 


(15) 


and  therefore  for  j  =  1,  2,  .  .  .  ,  mn,  we  have 


l°_1)(x)  =  £  ,  (x'l  i+1  f(x)ja>  +  ...  +  £  (x'“-i+1  f +  r®  (x)  .(16) 


(j> 


From  these  observations  and  equation  (9)  it  is  clear  that  as  x  and  n  become  large  G[F(x),  f(x),  .  .  .  , 

_  jt  \ 

'(x);  ajj(x)]  approaches  S  and  one  would  think  it  would  approach  S  more  rapidly  than  F(x). 
This  leads  us  to  the  following  definition,  which  is  essentially  just  an  economy  of  notation  over  our 
previous  remarks. 

Definition  L  Suppose  for  x  >  a,  that  i<”+-  *\x)  exists,  and  Hmn+1(cj  ,  a~(x))  ^  0,  then 

define  the  Gn  -transformation  by 


we 


Glm)[F(x);  a5J  =  g[f(x),  f(x),  ...J 


,(mn  — 1) 


(x);  aj  j(x) 


(17) 


The  G^^-transform  defined  by  (17)  is  in  fact  the  limiting  case  of  the  D-transformation  which  we 
spoke  of  in  the  introduction.  In  fact,  had  we  generated  the  system  of  equations  defined  by  equation 
(15)  by  incrementing  x,  i.e.,  using  F(x+kh),  instead  of  differentiating,  we  would  have  been  led  to  the 
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D-transformation  of  Levin  and  Sidi.  This  was  shown  in  [6].  To  obtain  Gn^  directly  as  a  limit  of  the 
D-transformation  one  begins  with  the  D-transform  and  takes  the  limit  as  h  — 0  after  several  row  and 

column  manipulations.  The  procedure  is  exactly  analogous  to  the  method  used  in  [2]  to  obtain  the 

(  ) 

Bm-transform  as  the  limit  of  the  Gm-transform.  It  follows  by  definition  that  Gj  =  Bm  . 


,(n») 


3.  Major  Results 

In  this  section  we  prove  rigorously  the  fundamental  properties  of  the  Gn  '-transformation 
suggested  by  our  previous  discussion.  The  first  theorem  needs  no  proof  and  is  really  just  a  restatement 
of  our  previous  remark  but  we  choose  to  state  it  formally  for  emphasis. 

Theorem  L  If  e(x)  satisfies  equation  (4)  with  •  =  a^j  =  .  .  .  =  am  j  =  0  when  i>n,  and 
Hmn+l(cj;  aij)  ^  °>  then  for  x  ^  a 


Gim)[F(x);  aj  J  =  j°°  f(u)du  . 


(18) 


,(m) 


The  value  of  a  transformation  such  as  Gq  does  not  lie  in  a  theorem  such  as  Theorem  1,  even  though 
such  theorems  are  of  interest.  However  it  is  the  next  theorem  which  justifies  our  interest  in  the  Gn  - 
transformation. 

(tu) 

Theorem  2.  Suppose  that  Gq  '|lF(x);  a- is  defined,  and  that  for  all  x  >  a,  f(x)  satisfies  equation  (4), 
and  that  m  >  0  is  the  smallest  integer  for  which  this  is  true.  Suppose  further  that  F(x)  — *  S  <  oo  as  x 
— ►  oo,  and  that  equation  (11)  holds.  Then 


S-G 


"m>[FW;  .  0  (,-») 

S-F(x)  UV  )■ 


(19) 
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Note  that  not  only  does  (19)  imply  that  Gn  [F(x);  a:  :]  converge  to  S,  but  that  it  does  so  more 

1 J 


rapidly  at  a  rate  proportional  to  x  n  as  x  — oo  . 


Proof.  Due  to  the  length  and  detail  required  to  prove  the  results  we  shall  divide  the  proof  into  two 
parts.  First,  we  will  show  the  results  in  some  detail  for  m=l.  Following  that  we  will  only  sketch  the 
proof  for  m  >  1. 

Case  I:  m  =  1. 

We  first  show  that  the  denominator  in  (8) 


tj  /  .  a  \  Anfn(x) 
Hn+l(cj’  aijJ  -  (n-l)n  ’ 


(20) 


I, 

for  some  An  ^  0.  Let  Pn(x)  =  Hn_^_2 (cj>  ajj)  an^  g(x)  =  x  f(x),  denoted  simply  by  g.  Then 


PnW  = 


a12(x) 

ai3(x) 

al,n+lW 

*22(x) 

^(x) 

a2,n+lW 

an,2<x) 

an,3(x) 

... 

an,n+l(x) 

g 


I 


kM!)'g'+(i)"g 


t  m 

j=oVjA 


(j)  (n-j) 
g 


,»-le"+2(xn1-1)  g'+(xn-l)  j|0('i)(x»-l 


0) 


(n-i) 


By  elementary  row  and  column  operations  we  obtain 


h 


1 


hn 


Pn(x)  =  gn~1 


where 

hj  =  g\ 

h2  =  g"-2g,h1/g  , 


g(»-j) 


hj/s- 


From  (16)  we  note  that  f(x)  -  a*j  x  T(x),  which  implies  that 
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and 


fG)(x)  -  dj  X  j£lf(x), 


(22) 


(23) 


for  some  constants  dj  and  dj,  j  =  1,  2,  ...  .  Note  that  dj  =  l/a*p  where  a|j  0  by  assumption  .  It 
then  follows  that 

(i-l)f, 

(24) 


hj  -  Cjf(x)/x  1  , 


for  some  constants  Cj,  j  =  1,  2,  ...,  n. 


(i).  If  £j  <  0,  then  it  is  easily  derived  from  (16)  that 

dj  =  dj  =  dj  ,  j  =  !,  2. . 

Then 


Cj  =  (-I)*  d|  ,  j  =  1,  .  .  .  ,  n 


Therefore 


pnW  -  g“  1  hn 


ay  ...  ay-11 


,  (n-i, 


n-n  •••Wn-1 


—  „n  —  1. 

=  g  hn 


(n  — l)n 


g 


n  — 1  cnf(x) 


bn  _  *  f°(x) 

(n  —  l)n  ‘  n  (n—  l)n 

X  X 


(25) 


where  An  =  cnbn  and  bn  is  a  constant.  It  is  seen  from  elementary  theory  of  differential  equations  that 
bn  ?£  0,  which  implies  that  An  ^  0. 

(ii).  Now  we  suppose  that  l  j  =  1.  Then  one  can  show  that 


10 


and  thus 


dj  =  (d1+l)  dj^-l)  .  .  .  (dj-j+2),  j  =  1,  2,  . 


•  •  ) 


=  An  g11-1^) 


(n-l)(n+l) 


(n-l)n  ’ 


where  An  =  bn  fj  (d,  — i+2)  for  some  bn  ^  0,  and  thus  An  ^  0  since  it  is  easily  shown  that 
i=l 

—  l<ajj<0.  We  have  therefore  proved  (20). 

We  now  consider  the  numerator  of  S  —  Gn  ^[F(x);  a;  •]  .  Denoting  the  numerator  by  Qn(x),  we 

lg 

have 


rn(x) 

rn(x) 

.S,n)(«) 

au(x) 

a12(x) 

al,n+l(x) 

QnW  = 

*2i(x) 

®22(x) 

^.n+lW 

an,l(x) 

an2(x) 

an.n+l(x) 

11 


where  hg  =  g,  Vq  =  rn(x)  and  v/s  correspond  to  hj’s  when  g  is  replaced  by  rn(x).  From  (14)  and  (22) 


we  know  that 


Thus,  by  (23), 

v0  =  rn(x)  =  0(x_n£(x)), 


=  rn(x)  =  O^x  “  1  e(x)J  , 

/  _n— 2 1,  \ 

v2  =  rn'(x)-2g'v1/g  =  Ol  x  1  e(x)  )  , 


We  then  have 


—  n  — n£ 
x 
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Qn(x)  =  'f0PnM-vl  « 


+  .  .  .  +  (-l)nvn  gn  1 


=  o(x  ”  <«)  P„W  +  O^V,  g"  \  ^n-DUn-l) 

+  8°“‘  h0  x(n-l)U„-a)  +  '  '  '  +  °(V"  *"'* 

=  o(x-nf(x))pn(x)  +  ofx  n+('  ^  1  ^(x)Pn(x) 


O^X  n+(n  2)(£l  1}  c(x)Pa(x)j  +  .  .  .  +  o(x-n  €(x)Pn(x)) 


=  C)(x  n«(x))pn(x)  . 


(26) 


Hence  by  (25)  and  (26),  we  obtain 


S-gS11)[F(x);  ay]  =  Qn(x)  /  Pn(x)  =  o(x“n£(x)) 


and 


S-Gi1)[F(x);  a}J]  _ 

- —  =  Ofx  nJ  ,  as  x— oo  , 


S-F(x) 


which  completes  the  proof  for  m=l. 


Case  II:  m  >  2. 


Because  the  proof  is  quite  lengthy  and  its  scheme  is  similar  to  Case  I,  we  only  sketch  it.  From 
(10),  we  can  asymptotically  express  /m\x)  in  terms  of  f(x),  f^x),  .  .  .  ,  /m  *\x)  as  follows: 


/m)«  -  £l  4  „  x  'l'0  /k)(x)  , 
k=l 


(27) 


where  q  are  nonzero  constants,  r^  q  =  £^  —  £m,  £q  =  0.  Because  the  relative  error  of  to  the 
right  hand  side  of  (27)  is  a  linear  combination  of  negative  powers,  we  have  the  following  lemma  whose 
proof  is  straightforward. 


Lemma  1.  For  q  =  1,  2,  ...  , 


k=0 
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(28) 


where  /?.  and  r.  are  constants,  and  r,  >  r.  ,  —  1. 

From  equation  (2),  /m+ci'l(x)  can  asymptotically  expressed  in  terms  of  f(x),  f^x), 

/m  ^(x).  Moreover,  one  can  show  that  there  exists  an  integer  ig  <  1  such  that  /q+*\x)  and 
xVq)(x)  have  the  same  order  when  x— ►co,  for  q  =  1,  2,  ...  .  To  be  specific,  suppose  m=2.  Let  (x) 
=  x  k  ^(x),  and  denote  it  by  g^,  k=l,  2  .  Then  and  x  ®gjp,  have  the  same  order,  for  j  = 

1,  2,  ...  .  Using  the  same  notation  as  Case  I,  we  have  the  denominator 


•22W  *23  (x) 


•2,2n+lW 


in,2(x)  an,3<X) 


an,2n+lW 


PnW  =  82 


an+2,2  an+2,3W 


n+2,2n-f  1' 


■2n,2n+lW 


There  are  (2n)!/(n!)  minors  of  order  n  obtained  from  the  first  n  rows.  Multiply  each  of  these  minors 
by  its  complementary  minor  and  by  a  proper  sign  factor.  The  sum  of  all  such  products  is  Pn(x)  by  the 
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Lapalace  expansion.  We  can  prove  that  each  product  is  at  most  of  order  f11  *(x)  (V(x)j 

—  n(n  —  l+(n+l)ig  — 

x  .In  fact,  ignoring  the  sign  factor,  the  first  such  product  is 


n+1 


ff(n) 

gl 


*2200  •  •  •  a2,n+l(x) 


i,2W 


an,n+l00 


(n+1) 

g2 


(2n) 

g2 


u+2,n+2 


(X)  .. 


11+2,211+1 


oo 


*2n,n+2W 


^n^n+lW 


b  (  g'1  1  " 

-  D0  ^  (n-l)i0  *1  (n-l)n/2; 


,  (n)  v 

r*2l  />>\“  1 

nin  lg2  )  (n  — l)n/2 
'x  U  x  ' 


-  bl  fQ-1(x)(f,(x))I1+l  .-“("-i-Kn+Dio-^-ej)  f 


where  bg  and  b^  are  constants,  by  employing  the  technique  used  in  the  proof  of  case  I.  Thus,  omitting 
the  detailed  tedious  arguments,  we  have 


PnOO 


Anf°  — X(x)  (^(x))D+1 


—  n(n—  l+(n+l)d  —  —  t^) 


where  An  is  a  nonzero  constant.  Similarly  to  the  proof  in  case  I,  we  can  also  show  that 


Qn(x)  =  0(x  n<r(x))  Pn(x). 


16 


The  case  of  m=2  is  then  proved.  The  proof  for  the  case  of  m>2  is  directly  analogous  to  the  case  m=2. 
It  is  however  extremely  laborious  and  we  do  not  include  it  here. 

Corollary  1.  Assume  the  general  conditions  in  Theorem  2.  If  max  I  aZ  n-i_i  I  ^  0,  then 

l<k<m 


S-ajffifrW; 

S-cir'ptx);  a;Jj 


n  =  1,  2, 


In  other  words 


r(m) 

Gn+1 


converges  to  S  faster  than 


at  a  rate  proportional  to  1/x. 


The  proof  is  straightforward  from  the  proof  of  Theorem  1. 


(29) 


,(m) 


4.  Applications 

In  Section  3  we  have  established  the  fast  convergence  rate  of  Gq  ''-transformation.  We  now 
consider  some  examples  to  show  the  numerical  accuracy  of  the  transformation. 


Example  l.  We  know  that 


oo 

r(a+l)  =  J  xa  e~x  dx  , 
0 


a  >  — 1  . 


The  integrand  f(x)  =  xae  x  satisfies 

f(x)  =  a^x  *V)  • 

Equation  (30)  implies  that  m=l  and  fj=0.  Rewriting  (30)  as 


(30) 


xf*(x)  =  (a-x)  f(x) 
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and  taking  the  r-th  derivative  on  both  sides  leads  to 


/r+1)(x)  =  (z=l-l^J  f(r)(x)  -  £  \)  . 


Using  (31),  the  Gn  -transform  can  be  carried  out  for  any  n.  To  be  specific,  let  a  =  |  .  Then  T(a+1) 
=  =  3.3233509704478  ...  .  Table  1  lists  the  comparison  of  the  G^- transform  and  Levin  and 

Sidi’s  Dq  ^-transform.  Adopting  their  notations  the  D^-transform  in  Table  1  uses  the  information  of 
F(x)  and  f(x)  at  x=l,  2,  ...  ,  n+1  while  the  G^  uses  F(x)  at  x  =  n-fl  and  the  derivatives  at  the  same 
point.  The  Gq  ^-transform  apparently  gives  better  approximations  in  this  case. 


Example  2.  In  the  second  example  we  follow  Example  4.4  in  [6].  The  goal  is  to  approximate  the 


integral 


oo 

/  dx  =  1.460362116753119547... 
J  1+x2 


The  integrand  f(x)  satisfies  the  differential  equation 


f  —  _  5x2+4x+l  f/  _  (x  +l)(x-H)  j, 
2(2x+l)  2(2x+l) 


Therefore  m=2,  i ^  =  1,  £«  =  2.  Rewriting  the  above  equation  as 


(4x+2)f(x)  =  —  (5x2+4x+ l^x)  -  (x^+x2+x+l)  f"(x) 


and  taking  r-th  derivative  on  both  sides,  we  obtain,  after  simplification, 


f<'+2)(x,=  -|[(3r+5)x2  +  (2r+4)  x+r+l]  /f+1)(x)  +  [(3r2+7r+4)  x+r2+3r+2]  f(r)(x) 


18 


+  r(r+l)2  *\x)j  /  (x3+x2+x+l) 


Thus  all  the  derivations  can  be  calculated  iteratively. 

(2)  (2)  (2) 

Table  2  compares  the  Gn  -transform  with  the  Dn  -transform,  where  Dn  use  F(x),  f(x)  and 

f*(x)  at  x  =  Xj  =  e^  ^  j=l,  •  •  .  ,  2n+l  and  uses  F(x),  f(x),  f^fx),  ...  ,  f^D+^\x)  at  x  = 

(2)  (2) 

x2n-t-l’  t^ie  lar8est  value  of  Xj’s.  Again  G^  is  far  better  than  Dq  .  Note  that  the  order  of  the 
determinants  involved  in  the  calculations  is  2n+l. 

An  important  application  of  the  G^^-transform  is  to  approximate  tail  probabilities.  In  this 
case,  calculations  of  the  transform  need  not  require  any  integration.  The  transformation  is  therefore 
useful  for  determining  approximation  functions  for  tail  probabilities.  Now  suppose  that  f(x)  is  a 
probability  density  function  (PDF)  and  we  wish  to  approximate  the  tail  probability 


oo 

P(a)  =  J  f(u)du  . 


/  \ 

If  we  let  a=x  in  (1)  then  F(x)  =  0,  so  that  it  is  clear  that  Gn  '  in  (17)  depends  only  on  f(x),  f^x),  ...  , 

imn+m  —  1).  .  .  .  ,  (m)r  ,, 

r  (x),  producing  an  approximation  function  for  P(a).  We  use  the  notation  Gn  |f(x)J  to 

indicate  this  specific  case.  This  method  has  been  shown  in  [5]  to  be  highly  accurate  even  with  small 

values  of  order  n  (n=l,  2,  3)  in  approximating  tail  probabilities.  In  such  lower  order  transformations, 

we  can  usually  express  the  resulting  approximations  in  very  simple  forms  so  that  the  calculations  of  the 

determinants  are  not  needed  and  a  hand  calculator  is  sufficient  to  compute  the  approximations. 


Example  J.  To  give  a  specific  example  we  now  consider  an  application  first  given  in  [5].  That  is, 
(m) 

G^  n  -transform  to  the  Pearson  family  which  include  such  well  known  PDF’s  as  the  normal,  \2, 
student  t,  F  and  others  (see  [5]).  For  the  Pearson  family  of  PDFs 
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«*)  =  b°  +-^t  b^  I'M  • 


Clearly,  m=l.  The  tail  probabilities  of  all  these  distributions  may  therefore  be  well  approximated  by 

G^[f(x);  ajj(x)].  For  those  distributions  with  finite  support,  the  method  can  still  be  employed  simply 

by  transformating  the  data  to  a  distribution  with  infinite  support. 

fm) 

For  larger  values  of  m  and  n  the  form  of  the  GV  n  -transform  given  by  equation  (17)  is  a 
convenient  one,  and  one  which  is  easy  to  calculate  on  a  micro-computer.  However  for  m=l  and  n  <  3 
the  transform  can  be  simplified  to  the  extent  that  it  is  easily  computed  on  a  hand  calculator. 
Specifically  we  have 


[f(x)] 


-X  f2(x) 
x  f'(x)  +  €jf(x) 


= _ x  f2(x)  A(x) _ 

2  U  U  x2B(x)-£1(£1-l)f2(x)-xf'(x)A(x)  ’ 


and 


G(lV(x)]  = 


xf2(x)  C(x) 


x3D(x)  +  3(£j  — 2)x2f(x) 


B(x)-(f'(x))2 


+  (£1-l)(«,-2)f2(x)E(x) 


where  A(x)  =  xf,(x)  +  2(£j  — l)f(x)  , 

B(x)  =  f(x)  f"(x)  -  (f'(x))2  , 

C(x)  =  3(£1-2)  f(x)[xf#(x)+(£1-l)f(x)]  -  x2|j2B(x)  -  (f^x))  , 

D(x)  =  -f2(x)  /3)(x)  +  6f*(x)  B(x)  , 

E(x)  =  -3xf*(x)  ~  «if(x)  • 
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To  demonstrate  the  precision  of  these  approximations,  we  consider  the  standard  normal  case. 
Finding  a  proper  functional  approximation  to  the  tail  of  a  normal  distribution  is  a  problem  which  has 
a  rich  history  in  statistics  and  many  approximations  have  been  proposed.  In  some  applications  one 
requires  a  tail  probability  that  is  accurate  in  the  extreme  tails  as  well  as  at  the  usual  nominal 
significance  levels.  This  problems  has  recently  been  addressed  by  Hawkes  [7]  who  employed  an  ad  hoc 
approach  to  obtain  a  specialized  approximation  which  is  accurate  in  the  very  extreme  tails.  If  f  is  a 
standard  normal  density  function,  i.e., 


,,  ,  1  —  x2/2 

f(x)  -  e 


then 


and  clearly  i ^  =  —1.  Furthermore 


f(x)  =  -  i  f'(x) 


f^M  =  -J1  1)(x)  +  (l-k)f(r  2)(x),  r  >  2, 


,(1) 


so  that  Gn  [f(x)]  is  easy  to  calculate  for  any  n.  Moreover  it  is  easily  derived  from  the  definition  that 


.(Dr 


and 


Gi  '[f<x)]  =  -f-  f(x) 
x  +1 


G^Wx)]  =  -T-~54)- - fW 

(xj+1)(x2+4)-2 


G^[f(x)]  =  -^-f+2)(f+9)  f(x) 

x2(x2+3)(x"+9)+6 


Table  3  compares  G^  \  G^  ^  and  Gg  ^  with  the  approximations  and  given  in  [7];  see  [7]  for 
details.  Note  that  G^Vw]  is  essentially  as  good  as  the  best  of  and  Qjj0  and  G^^x)]  is  better. 
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In  fact,  these  new  approximations  are  exact  to  many  more  significant  digits  than  those  shown  in  Table 
3  in  the  extreme  tails.  For  example,  when  x=8,  the  absolute  eriors  of  Go  ^  and  Gj  ^  are  3.4  X  10~ 
and  2.5  X  10~“^  respectively  while  those  of  and  Q{J2  are  ^  ^  and  2.5  X  10 

respectively.  Notice  that  all  G^[f(x)],  G^[f(x)]  and  G^[f(x)]  are  a  rational  function  times  the 
density  function.  This  is  true  for  all  other  distributions  considered  in  [5]. 


5.  Concluding  Remarks 

In  this  paper  we  have  extended  the  Bn-transformation  in  manner  exactly  analogous  to  the  Levin- 
Sidi  generalization  of  the  Gn-transformation.  As  in  the  case  of  the  Levin-Sidi  lesults,  the  consequences 
of  such  an  extension  are  substantial,  <»s  was  demonstrated  by  Theorem  2  and  the  examples.  It  should 
be  noted  that  the  examples  given  are  inonotonic  for  large  x.  This  is  certainly  not  necessary  for  the 
G^^-transform  to  be  effective.  However  as  might  be  expected  from  our  motivation  of  the  G^^- 
transform,  these  are  the  types  of  functions  for  which  the  G^^-transform  would  be  expected  to  perform 
the  best.  For  example,  the  G^™^ -transform  can  be  used  to  accelerate  the  convergence  of  the  integral 


oo  . 


n(2) 


J  5iB_t  dt  but  is  not  as  effective  as  Levin  and  Sidi’s  D'n-transfoim. 
0 
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Table  1 


Errors  of  the  Approximations 


and  D 


(1) 


to 


dx  =  3.32335097044784227.  .  .  . 


n 

|  G^W+l);  a-J  -  r(7/2)  | 

1  E>n1)[F(l);  1]  -  1 

4 

1.7X10- 3 

5.7X10-2 

6 

1.7X10- 6 

JO 

Oi 

X 

1 — • 
o 

1 

8 

3.1X10-9 

2.4X10“  6 

10 

7.3X10“ 12 

2.8X10-8 

12 

1.9X10-14 

3.6X10-10 

14 

2.9X10-16 

5.1X10"12 
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Table  2 


(2)  (2) 

Errors  of  the  Approximations  Gn  and  Dn  to 


oo 

-I 


log  X 
1+x2 


dx  =  1.460362116753119547.  .  .  . 


n 

1  GpIFtW  *ijl  -  S  1 

1  D<2)[F(x.)]  - 

2 

8.1X10- 3 

3.2X10“ 1 

3 

7.6X10“ 4 

8.2X10“3 

4 

3.6X10-5 

5.0X10-4 

5 

5.0X10-7 

3.9X10“4 

6 

3.6X10-11 

6.3X10“5 

7 

1.7X10“ 12 

2.9X10“6 

8 

1.1X10”15 

3.4X10-8 
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Table  3.  Approximations  to  the  upper  taii  of  the  normal  distribution. 


X 

true  tail 

QL2 

QH2 

o',1’ 

«41’ 

c(1) 

G3 

1.2 

.11507 

.11737 

.11402 

.09550 

.11244 

.11504 

1.6 

.05480 

.05589 

.05465 

.04985 

.06452 

.05486 

2.0 

.02275 

.02314 

.02274 

.02160 

.02273 

.02276 

2.4 

.008198 

.008312 

.008197 

.007951 

.008199 

.008200 

3.0 

.001350 

.001364 

.001350 

.001330 

.001350 

.001350 

4.5 

3.398x10“® 

3.414x10“® 

3.396x10“ 6 

3.385x10“® 

3.398x10“® 

3.398x10 

6.0 

9.866x10“*® 

9.891x10“*® 

9.861xl0“10 

9.853x10“ 10 

9.866x10“ 10 

9.866x10 

8.0 

6.221x10“ 16 

6.229x10“*® 

6.218x10“ 16 

6.218xl0“16 

6.221x10“*® 

6.221x10 

10.0 

7.620xl0“24 

7.625xl0~24 

7.618xl0“24 

7.618xl0“24 

7.620xl0“24 

7.620x10 

14.0 

7.794xl0“45 

7.796xl0~45 

7.793x10  “45 

7.793x10  “45 

7.794x10  ”45 

7.794x10 

18.0 

9.741xl0~73 

9.742x10“  73 

9.741xl0“73 

9.741xl0“73 

9.741xl0“73 

9.741x10 
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